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Abstract 

We study classical Hamiltonian systems in which the intrinsic proper time evolution pa- 
rameter is related through a probability distribution to the physical time, which is assumed 
to be discrete. 

This is motivated by the "timeless" reparametrization invariant model of a relativistic 
particle with two compactified extradimensions. In this example, discrete physical time 
is constructed based on quasi- local observables. 

Generally, employing the path-integral formulation of classical mechanics developed by 
Gozzi et al., we show that these deterministic classical systems can be naturally described 
as unitary quantum mechanical models. The emergent quantum Hamiltonian is derived 
from the underlying classical one. It is closely related to the Liouville operator. We 
demonstrate in several examples the necessity of regularization, in order to arrive at 
quantum models with bounded spectrum and stable groundstate. 



*Lecture given at DICE 2002. To be published in: Decoherence and Entropy in Complex Systems, 
Lecture Notes in Physics (Springer- Verlag, Berlin 2003). 
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1 Introduction 

Recently we have shown that for a particle with time-reparametrization invariant dynamics, be 
it relativistic or nonrelativistic, one can define quasi-local observables which characterize the 
evolution in a gauge invariant way |H |2] . 

We insist on quasi-local measurements in describing the evolution, which respect reparametriza- 
tion invariance of the system. Then, as we have argued, the physical time necessarily becomes 
discrete, its construction being based on a Poincare section which reflects ergodic dynamics, 
by assumption. Most interestingly, due to inaccessability of globally complete information on 
trajectories, the evolution of remaining degrees of freedom appears as in a quantum mechanical 
model when described in relation to the discrete physical time. 

While we pointed out in the explicitly ergodic examples of references P3 12] that such emer- 
gent discrete time leads to what may, for obvious reasons, be called "stroboscopic" quantization, 
we report here how this occurs quite generally in classical Hamiltonian systems, if time is dis- 
crete and related to the proper time of the equations of motion in a statistical way (3J. In 
our concluding section, we will briefly comment about extensions, where the prescribed prob- 
abilistic mapping of physical onto proper time shall be abandoned in favour of a selfconsistent 
treatment. A closed system has to include its own "clock", if it is not entirely static, reflecting 
the experience of an observer in the Universe. 

Previous related work on the "problem of time" has always assumed that global features 
of the trajectory of the system are accessible to the observer. This makes it possible, in 
principle, to express the evolution of an arbitrarily selected degree of freedom relationally in 
terms of others jU |3] . Thereby the Hamiltonian and possibly additional constraints have been 
eliminated in favour of Rovelli's "evolving constants of motion" 6J. For a recent development 
aiming at solving the constraints after discretization see [7j. While appealing by its conceptual 
clarity, incorporating nonlocal observations seems unrealistic to us in any case. 

In distinction, we point out that the emergent discrete time in our approach naturally leads 
to the "stroboscopic" quantization of the system [H 121 13] • Quantum theory thus appears to 
originate from "timeless" classical dynamics, due to the lack of globally complete information 
0- 

Another approach to deterministically induced quantization is proposed in [H], where the 
consequences of incomplete statistics are analyzed, leading towards Euclidean quantum field 
theory under very general assumptions. Various other arguments considering quantization as an 
emergent property of classical systems have recently been proposed, for example, in references 
[Jill HH U21 EE] j concerning quantum gravity and dissipation at a fundamental level, "chaotic 
quantization" , and matrix models, respectively. 

Our approach tries to illuminate from a different angle how to arrive at quantum models 
which describe dynamically evolving systems. In particular, we believe that there may be an 
intimate connection with how the "problem of time" is resolved for a local observer, namely by 
counting suitably defined and locally measurable incidents. 

We remark that the possibility of a fundamentally discrete time (and possibly other discrete 
coordinates) has been explored before, ranging from an early realization of Lorentz symmetry 
in such a case |T3] to detailed explorations of its consequences and consistency in classical 
mechanics, quantum field theory, and general relativity (T^l E3 El- However, no detailed 
models giving rise to such discreteness have been proposed. Quantization, then, is always 
performed in an additional step, as usual. 
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In particular, the work by Gambini, Pullin et al. aims at a consistent canonical quantiza- 
tion of gravity via discretization [7JEI- Discretization of time is performed in a static fashion, 
i.e. independently of the evolution. As shown there, the major advance lies in the possibility 
to satisfy the constraints, in principle, by suitably choosing the Lagrange multipliers. How- 
ever, the extraneous discretization is reflected by the persistence of a discrete variable n after 
quantization, which apparently has no physical meaning. 

Presently, we shall see a vague resemblance to the persistence of proper time r in our 
approach, as long as the discrete physical time is related through a given probability distribution 
to proper time. Our formalism, however, is set up in such a way that the "clock" degrees of 
freedom can be treated dynamically as part of the system. This should help to pinpoint the 
role of proper time in the resulting model where, in our case, quantization emerges instead of 
being imposed on the system. 

In the following section, we recall the model from [2], in order to motivate the emergence of 
discrete time from quasi-local observations. This is the starting point of our heuristic derivation 
of a quantum mechanical picture of what appear fundamentally classical systems. 

To put our approach into perspective, we remark that there is clearly no need to follow such 
construction leading to a discrete physical time in ordinary mechanical systems or field theories, 
where time is an external classical parameter, commonly called "t". However, assuming for 
the time being that truly fundamental theories will turn out to be diffeomorphism invariant, 
adding further the requirement of the observables to be quasi-local (modulo a fundamental 
length scale), when describing the evolution, then such an approach seems natural, which may 
lead to quantum mechanics as an emergent description or "effective theory" on the way. 

2 Discrete Time of a Relativistic Particle 
with Extradimensions 

We consider the (5+l)-dimensional model of a "timeless" relativistic particle (rest mass m) 
with the action: 



Here A stands for an arbitrary "lapse" function of the evolution parameter s, x M = dx^/ds (/i = 
0, 1, . . . , 5), and the metric is g^ v = diag(l, —1, . . . , —1). Units are such that h = c = 1. 

With this form of the Lagrangian, instead of the frequently encountered Lds oc (g^dx^dx") 1 ! 2 
which emphasizes the geometric (path length) character of the action, the presence of a con- 
straint is immediately obvious, since there is no s-derivative of A. 

Two spatial coordinates, x 4,5 in (0), are toroidally compactified: 




(1) 



where the Lagrangian is defined by 




(2) 



x 



,4,5 = 



2nR[<f> 4 > 5 } , 



(3) 
(4) 




(j) e [n, n+ 1[ , 



for any integer n, i.e. the angular variables are periodically continued; henceforth we set R = 1, 
for convenience. Alternatively, we can normalize the angular variables to the square [0, l[x [0, 1[, 
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of which the opposite boundaries are identified, thus describing the surface of a torus with main 
radii equal to one. 

While full Poincare invariance is broken, as in other currently investigated models with 
compactified higher dimensions, the usual one remains in fourdimensional Minkowski space 
together with discrete rotational invariance in the presently two extradimensions; also trans- 
lational symmetry persists. Furthermore, the internal motion on the torus is ergodic with an 
uniform asymptotic density for almost all initial conditions, in particular if the ratio of the 
corresponding initial momenta is an irrational number. 

Setting the variations of the action to zero, we obtain: 

^ = \\~ 2 x^-m 2 ) =0 , (5) 



5X 2 
5S d 



6x^ ds 



— (A x M ) = . (6) 



In terms of the canonical momenta, 

the equations of motion © become simply p^ = 0, while (0) turns into the mass-shell constraint 
p 2 — m 2 = 0. 

The equations of motion are solved by: 

x »(s) =x?-p» f ds'A(s') = xf + P »t(s) , (8) 
Jo 

where the conserved (initial) momentum p M is constrained to be on-shell and Xi denotes the 
initial position. Here we also defined the fictitious proper time (function) r, which allows 
us to formally eliminate the lapse function A from (jHJ)-©, using x^{t) = x IJ, {s) and i^(s) = 

-\(s)d T X»(T). 

In order to arrive at a physical space-time description of the motion, the proper time needs 
to be determined in terms of observables. In the simplest case, the result should be given by 
functions x M ^ (x°), provided there is a physical clock measuring x° = X® +p°r. 

Similarly as in the nonrelativistic example studied in [Q. the lapse function introduces 
a gauge degree of freedom into the dynamics, which is related to the reparametrization of 
the evolution parameter s. In fact, the action, (Q)-©, is invariant under the set of gauge 
transformations: 

s = f(s') , af(s) = s"V) , AOO^AV) ■ (9) 

It can be shown that the corresponding infinitesimal tranformations actually generate the evo- 
lution of the system. This is the basis of statements that there is no time in systems where 
dynamics is pure gauge, i.e. of the "problem of time". We refer to PP for further discussion. 

Instead, with an evolution obviously taking place in such systems, we conclude from these 
remarks that the space-time description of motion requires a gauge invariant construction of a 
suitable time, replacing the fictitious proper time r. To this we add the important requirement 
that such construction should be based on quasi-local measurements, since global information 
(such as invariant path length) is generally not accessible to an observer in more realistic, 
typically nonlinear or higherdimensional theories. 
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2.1 "Timing" Through an Extradimensional Window 

Our construction of a physical time is based on the assumption that an observer in (3+1)- 
dimensional Minkowski space can perform measurements on full (5+l)-dimensional trajectories, 
however, only within a quasi-local window to the two extradimensions. In particular, the 
observer records the incidents ("units of change") when the full trajectory hits an idealized 
detector which covers a small convex area element on the torus (compactified coordinates 
x 4 - 5 ). 1 

Thus, our aim is to construct time as an emergent quantity related to the increasing number 
of incidents measured by the reparametrization invariant incident number: 

1 = H ds'X(s')D(x 4 (s'),x 5 (s')) , (10) 

J Si 

where x 4,5 describe the trajectory of the particle in the extradimensions, the integral is taken 
over the interval which corresponds to a given invariant path and the function D 

represents the detector features. Operationally it is not necessary to know the invariant path, 
in order to count the incidents. 

In the following examples we choose for D the characteristic function of a small square of 
area d 2 , D(x 4 ,x 5 ) = Cd(x 4 )Cd{x 5 ), with Cd(x) = Q(x)(l — Q(x — d)), which could be placed 
arbitrarily. Our results will not depend on the detailed shape of this idealized detector, if it is 
sufficiently small. More precisely, an incident is recorded only when, for example, the trajectory 
either leaves or enters the detector, or according to some other analogous restriction which could 
be incorporated into the definition of D. Furthermore, in order not to undo records, we have to 
restrict the lapse function A to be (strictly) positive, which also avoids trajectories which trace 
themselves backwards (or stall). The records correspond to a uniquely ordered series of events 
in Minkowski space, which are counted, and only their increasing total number is recorded, 
which is the Lorentz invariant incident number. 

Considering particularly the free motion on the torus, solution (jHJ) yields: 

[0(r)] = [0 o + 7rr] , (11) 

where (f) is the vector formed of the angles 4 ' 5 , and correspondingly fr, with ir 4,5 = p ' 5 /27rR; the 
quantities in (fTTj) are periodically continued, as before, see (EH)-®- Without loss of generality 
we choose O = and 7r 5 > 7r 4 > 0, and place the detector next to the origin with edges aligned 
to the positive coordinate axes for simplicity. 

Since here we are not interested in what happens between the incidents, we reduce the 
description of the internal motion to coupled maps. For proper time intervals At with 7r 4 ■ At = 
1, the </> 4 -motion is replaced by the map m — > m + 1, where m is a nonnegative integer, while: 

[0 5 ] = [Pm] , (12) 

with P = 7r 5 /7r 4 > 1. Then, also the detector response counting incidents can be represented 
as a map: 

J(m + l) = /(m) + e(5-[0 5 ]) + e([0 5 ]-(l-^)) , (13) 

^^One could invoke a popular distinction between brane and bulk matter as in string theory inspired high- 
crdimcnsional cosmology, in order to construct more realistic models involving local interactions. 
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with 1(0) = 1, and where 5 = d/2nR corresponds to the detector edge length d, assumed to 
be sufficiently small, PS <C 1. The two G-function contributions account for the two different 
edges through which the trajectory can enter the detector in the present configuration. 

The nonlinear two-parameter map (|13|) has surpring universal features, some of which we 
explored in |2j. Here, first of all, following the reparametrization invariant construction up to 
this point, we identify the physical time T in terms of the incident number / from (fH3|) : 



T = 



5(7T 4 + 7T 5 ) 



(14) 



A statistical argument for the scaling factor 5 1 (tc a + tt 5 ) l , based on ergodicity, has been given 
in pQ, which applies here similarly 



o 
o 
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Figure 1: The physical time T as a function of proper time r with detector parameter 5 = .005 
and ratio of initial internal momenta P = a/31 (top), e, a/2, 7r (bottom) (see main text); upper 
two curves displaced upwards by +10 and +20 units, respectively, for better visibility (from 

ED- 



We show in Fig. 1 how the physical time T typically is correlated with the fictitious proper 
time t. The proper time is extracted as those m- values when incidents happen: r = m + 1 <^=^ 
I(m + 1) — I(m) = 1, corresponding to a particularly simple specification of the detector 
response. For a sufficiently small detector other such specifications yield the same results 
as described here. This is achieved by always rounding the extracted proper time values to 
integers, involving negligible errors of order 5,5/P <^ 1. 

We find that the time T does not run smoothly. This is due to the coarse-grained description 
of the internal motion: as if we were reading an analog clock under a stroboscopic light. In 
our construction, it is caused by the reduction of the full motion to a map (Poincare section), 
corresponding to the recording of the physical incidents by the quasi-local detector. 
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Furthermore, already after a short while, i.e. at low incident numbers, the constructed time 
approximates well the proper time r on average. 2 The fluctuations on top of the observed linear 
dependence result in the discreteness of the constructed time. 

While in reference [2] we further analyzed the statistical properties of the map considered in 
this example, we take from here only the result that a physical time can be constructed based 
on suitable localized observations. Furthermore, after embarking on some useful formal devel- 
opments in the next section, we will incorporate the resulting probabilistic mapping between 
discrete physical time and proper time of the equations of motion in Section 4. 



3 Classical Mechanics in Path-integral Form 

Classical mechanics can be cast into path-integral form, as originally developed by Gozzi, Reuter 
and Thacker ^H] , and with recent addenda reported in • While the original motivation has 
been to provide a better understanding of geometrical aspects of quantization, we presently use 
it as a convenient tool. We refer the interested reader to the cited references for details, on the 
originally resulting extended (BRST type) symmetry in particular. We suitably incorporate 
time-reparametrization invariance, assuming equations of motion written in terms of proper 
time. 

Let us begin with a (2n)-dimensional classical phase space Ai with coordinates denoted 
collectively by <p a = (q 1 , . . . , q n ; p 1 , . . . ,p n ), a = 1, . . . , 2n, where q,p stand for the usual 
coordinates and conjugate momenta. Given the proper-time independent Hamiltonian H(tp), 
the equations of motion are: 

^ = ^A. HM , (15) 

where oj ab is the standard symplectic matrix and r denotes the proper time; summation over 
indices appearing twice is understood. 

To the equation of motion we add the (weak) Hamiltonian constraint, Ch = H(ip) — e ~ 0, 
with e a suitably chosen parameter. This constraint has to be satisfied by the solutions of the 
equations of motion. Generally, it arises in reparametrization invariant models, similarly as 
the mass-shell constraint in the case of the relativistic particle fl . It is necessary when the 
Lagrangian time parameter is replaced by the proper time in the equations of motion. In this 
way, an arbitrary "lapse function" is eliminated, which otherwise acts as a Lagrange multiplier 
for this constraint. 

We remark that field theories can be treated analogously, considering indices a, b, etc. as 
continuous variables. 

Starting point for our following considerations is the classical generating functional, 

Z[J] = J H Vcp 5[ip a (r) - ^(r)] exp(* J dr J aV a ) , (16) 

where J = { J a =i,...,2n} is an arbitrary external source, S[-] denotes a Dirac ^-functional, and <p c i 
stands for a solution of the classical equations of motion satisfying the Hamiltonian constraint; 
its presence is indicated by the subscript "if" on the functional integral. The relevant boundary 
conditions shall be discussed in the following section. It is important to realize that Z[0] gives 



2 The apparent excursion for the parameter value P = ir in Fig. 1, does not persist for longer times, as shown 
m 12 
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weight 1 to a classical path satisfying the constraint and zero otherwise, integrating over all 
initial conditions. 

Using the functional equivalent of S(f(x)) = \df /dx\~^ ■ S(x — x ) , the 5-functional under 
the integral for Z can be replaced according to: 

5{ V a (r) - p« (r)] => 5[d T (p a - cu ab d b H] det[5^d T - u ac d c d b H] , (17) 

slightly simplifying the notation, e.g. d b = d/d(p b . Here the modulus of the functional deter- 
minant has been dropped ^Sl EH] • 

Finally, the 5-functionals and determinant are exponentiated, using the functional Fourier 
representation and ghost variables, respectively. Thus, we obtain the generating functional in 
the convenient form: 

Z\J] = J^VtpVWcVc exp (i J dr(L + J a ^ a )) , (18) 

which we abbreviate as Z[J] = f H T>& exp(i / drLj). The enlarged phase space is (8n)- 
dimensional, consisting of points described by the coordinates (y? a , A a , c a , c a ). The effective 
Lagrangian is now given by jlH| ITH] : 

L = \ a (d T ip a - uj ab d b H) + ic a (5 a b d T - u ac d c d b H)c b , (19) 

where c a , c a are anticommuting Grassmann variables. We remark that an entirely bosonic 
version of the path-integral exists ^J] ■ 

This completes our brief review of how to put (reparametrization invariant) classical me- 
chanics into path-integral form. 

4 From Discrete Time to "States" 

We recall from our previous example that the discrete physical time t has been obtained by 
counting suitably defined incidents, i.e., coincidences of points of the trajectory of the system 
with appropriate detectors |UI2]- Thus, it is given by a nonnegative integer multiple of some 
unit time, t = nT. Then, we would like to express the proper time r which parametrizes the 
evolution in terms of t. 

Here we assume instead that the physical time t is mapped onto a normalized probability 
distribution P of proper time values r: 

P(r; t) = exp ( - S(r; t)) , j dr P(r; t) = 1 . (20) 

For uniqueness, we require that if S(r;ti) and S , (r;t 2 ), for t\ ^ t 2 , have overlapping support, 
then they should coincide in this region. 

Thus, we describe the idealized case that the system can be separated into degrees of 
freedom which are employed in the construction of a physical "clock" , yielding the values of t, 
and remaining degrees of freedom evolving in proper time. Neglecting the interaction between 
both components, and the details of the clock in particular, we describe the relation between 
physical and proper time by a probability distribution. In this situation, the Hamiltonian 
constraint only applies to the remaining degrees of freedom, while generally the system will be 
constrained as a whole. These aspects were exemplified in detail in the simple models of PUI2]- 
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Correspondingly, we introduce the modified generating functional: 

Z[J] = J h dndTf J V$ exp (i jP dr L 7 - S^U) - S(r f ;t f )) , (21) 

instead of (fT%|) . using the condensed notation introduced there. In the present case, Z[0] sums 
over all classical paths satisfying the constraint with weight P^;^) • P(ry;ty), depending on 
their initial and final proper times, while all other paths get weight zero. In this way, the 
distributions of proper time values Tij associated with the initial and final physical times, tj 
and tf, respectively, are incorporated. 

Next, we insert 1 = / drP(r; t) into the expression for Z, with an arbitrarily chosen physical 
time t > t , and with ti = tf = t . We require the two sets of trajectories created in this way 
to present branches of forward (">") and backward ("<") motion. This leads us to factorize 
the path-integral into two connected ones: 

Z[J] = J dr P(r; t) ■ j dr f J exp (i j*' dr' L< - S(r f ; t )) 

■ j^j^ X>$> exp (i jT dr" L> - S(r 4 ; t j) 

■ II%>W-^)^W) ' ( 22 ) 

a 

where e(a < n) = 1, e(a > n + 1) = —1, and J = J > ,J < , depending on the branch. 
The ordinary ^-functions assure continuity of the classical paths in terms of the coordinates 
q a , a — 1, ... ,n, and reflect the momenta p a , a = 1, . . . , n, at proper time r. 

We observe that the generating functional will only be independent of the physical time t, 
in the absence of an external source, if we assume that the probability distribution P does not 
explicitly depend on time, — log P(r;t) = S(r;t) = S(r — t), and if we suitably specify the 
boundary conditions. We set: 

^ a > (r i ) = e(a)^(r / ) = a (t o ) , a = l,...2n . (23) 

Note that the boundary conditions are defined at the physical time to, to which correspond 
the distributed values of the proper times r^j. This establishes a one-to-one correspondence 
between both sets of trajectories. They could be viewed as closed loops with reflecting boundary 
conditions at both ends, to and t, and fixed initial condition at to- 
Exponentiating the 5-functions via Fourier transformation, the generating functional can 
be recognized indeed as a scalar product of a "state" and its adjoint. We define the normalized 
states by the path-integral: 

|r,vr a ;t) EEZfJ]- 1 / 2 ! dr^ exp (z £ + * dr' Lj - S(r,; t ) + m aV a {r + t)) , (24) 

and, similarly, the adjoint states: 

{r^ a -t\=Z[J}-^ 2 j dr f J h V§ exp (i J^dr' Lj - S(r f ;t ) - m a cp a (T + t)), (25) 

where the paths are forward and backward going as indicated by the integral boundaries in the 
exponent respectively, dropping ">, <"; note that the summation is to be read as J2a e ( a ) 7l af a 
in (EH). 
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The redundancy in designating the states, which depend on the sum of proper and physical 
time, only arises here, since the probability distribution P is assumed not to be explicitly 
depending on the physical time, for simplicity. 

The scalar product of two such states is now defined, and calculated, as follows: 

<*2|*i> = / drdTT P(t)<t, 7r;t a |r,7r;ti) = 5 t2jh , (26) 

with dir = TiaidTT a /27r). In particular, we have (t\t) = 1, which corresponds to using 
definitions (|24|) and (|25j) . Furthermore, for t\ ^ t 2 ? we find that states are orthogonal, (i 2 = 
0, by the symmetry of the motion on the forward and backward branches, for a correspondingly 
symmetric source J, and by uniqueness of the Hamiltonian flow generating the paths. 

A remark is in order here concerning the integration over d7r above, which originates from 
exponentiating the ^-functions of (|22|1 . Due to the presence of the Hamiltonian constraints on 
both branches of a trajectory, one of the ^-functions is redundant. We absorb the resulting 5(0) 
in the normalization of the states. 

Finally, the symmetry between the states and the adjoint states, given the stated assump- 
tions, is perfect. We find: 

(T,7r;t| = |<r,7r;t>* , (27) 

which is a desirable property of states in a Hilbert space (in "r, 7r-representation"). However, 
from our heuristic discussion this appears as a restriction which could be relaxed, resulting in 
a less familiar relation between the vector space and its dual. 3 



5 Unitary Evolution 

Following the same approach which led to the definition of states in ()24|) . (|25 p . we consider the 
time evolution of states in the absence of a source, J = 0. Suitably inserting "1", as before, 
and splitting the path-integral, we obtain: 

|r', tt'; f) = J drdn P(t)U(t\ tt'; t'\r, tt; t)\r, vr; t) , (28) 

with the kernel: 

/rr'+t' 
X?$ exp (i J dr" L + m' ■ ip(r' + *') - iir ■ ip(r + t)j , (29) 

where the integral is over all paths running between r + 1 and r' + 1', subject to the constraint; 
here we abbreviate Tt-ip = 7r a (p a . We interpret this as a matrix element of the evolution operator 
U(t'\t). 

Then, it is straightforward to establish the following composition rule: 

U(t"\t') ■ U(t'\t) = U(t"\t) , (30) 

where integration over the intermediate variables, say t',tt' , with appropriate weight factor 
P(r'), is understood. These integrations effectively remove the "1", which is inserted when 

3 One might consider only forward going paths, for example, with boundary conditions on the coordinates 
set at ±to- In this case, however, it is not obvious how to obtain the correspondent of l|27|) . 
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factorizing path-integrals, and link the endpoint coordinates of one classical path to the initial 
of another. 

Since the Hamiltonian constraint is a constant of motion, there is no need to constrain 
the path integral representing the evolution operator. Integrating over intermediate variables 
removes all contributions violating the Hamiltonian constraint, provided we work with properly 
constrained states. This will be further discussed in the following section. 

The physical-time dependence of the evolution operator amounts to translations of proper 
time variables. Therefore, we may study further properties of U without explicitly keeping 
it. This simplicity, of course, is related to the analogous property of the states, which we 
mentioned. 

We begin by rewriting the functional integral of (|29|) : 

U{r', tt'; t, tt) = J Vp 6[p a (r) - ip a d {r)} exp (V ■ p{r') - m ■ ip{r)) , (31) 

cf. Section 3, where the paths run between r and r', integrating over all initial conditions. 
Fixing the initial condition of a classical path, we can pull the exponential factors out of the 
integral, due to the 5-functional, and integrate over all initial conditions in the end: 

U(r' 1 n / ;r,n) = j dip^r) exp (in' • ip f {r') - in ■ ipi{r)) J Vp 6[p a (r) - p a cl (r)}, (32) 

where Pf{r') denotes the endpoint of the path singled out by the particular initial condition, 
ipi(r). The functional integral equals one. Then, we obtain the simple but central result: 

U (t', 7r'; r, 7r) = / dp exp (in exp[£(r' — r)] • ip — in ■ ipj (33) 
= Stf^-y-T) , (34) 
where dip = riad<^ a , an d with the Liouville operator: 

£ = -^-^' (35) 

which is employed in order to propagate the classical solution from the initial condition at r to 
proper time r'; uj is the symplectic matrix. 

Using one readily confirms (J3T?|) once again. In particular, then U(t\t')-U(t'\t) = U(t\t), 
which is not diagonal, in general, in this r, 7r-representation. We have: U(r', tt'; t\r, tt; t) = 
S(n', tt; t' — t), as defined in 

In order to proceed, we consider the time dependence of the evolution kernel £. It is 
determined by the equation: 

id T £(n', 7r; r) = — J dip exp (in' ■ ip(r) — n ■ <p^ tt' ■ u ■ —H{<p{r)) 

= H{n\-id^)£{n\n;T) , (36) 

with the effective Hamilton operator. 

— d 

H(n, -id n ) = tt ■ uj —H(ip)\ v= _ idiT . (37) 
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Here we used together with the equation of motion (JTSJ). The initial condition is: 

E(tt', tt; 0) = (27r) 2 "5 2n (7r / - tt) , (38) 

as read off from (J33J). 

Using (JHHJ), we finally obtain the Schrddinger equation which describes the evolution of the 
states in physical time: 

id t {T,7r\V(t)) = J dr'dn' P(T')id t £(TT,ir';T + t - t')(t' ,ir'\^(0)) 

= W(7r,-i9 ff )(r,7r|*(t)) . (39) 

Here we implicitly employed the relation (r, ^{^(t)) = (r + 1, 7r|\I/(0)} = (0, 7r\^(t + r)), in 
order to analytically continue to real values of t and to perform the derivative, despite that the 
physical time is discrete. 

Clearly, the hermitian Hamiltonian must be incorporated in a unitary transfer matrix, in 
order to describe the evolution through one discrete physical time step. It is plausible that 
presently the need for regularization of the Hamilton operator, demonstrated in subsequent 
sections, arises here. It is also conceivable that in a more realistic situation, with clock degrees 
of freedom forming dynamically part of the system, this particular complication is alleviated. 

Furthermore, considering stationary states, we have: 

(r,7r|* B (*)) = exp(-iM)<r,7r|*(0)) = exp(-iE(t + r))(0,7r|*(0)) (40) 
= exp(-i£(t + T))(7r|* E > , (41) 

due to the previously discussed additivity of proper and physical time in the present context. 
Similarly, the Hamiltonian TC is independent of the probability distribution P, mapping physical 
to proper time, since in the presently idealized situation the clock is decoupled from the system. 

Note that there is no H in our equations. If introduced, it would merely act as a conversion 
factor of units. On the other hand, there is an intrinsic scale corresponding to the clock's 
unit time interval T, which could be analyzed in a more complete treatment where clock and 
mechanical system are part of the Universe and interact. 

Before we will illustrate in some examples the type of quantum Hamiltonians that one 
obtains, we have to first address the classical observables and their place in the emergent 
quantum theory, in particular we need to implement the classical Hamiltonian constraint. We 
recall that in a reparametrization invariant classical theory the Hamiltonian constraint is an 
essential ingredient related to the gauge symmetry one is dealing with. 



6 Observables 

It follows from our introduction of states in Section 4, see particularly (|2i p -([25| ) . how the clas- 
sical observables of the underlying mechanical system can be determined. Considering observ- 
ables which are function(al)s of the phase space variables (p, the definition of their expectation 
value at physical time t is obvious: 

(0[<p];t) = JdTP(r;t)O[-i-J-r}logZ[J}\j =0 (42) 
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: J drdn P(r -t)(T,n;0\O[ip(r)]\r,n;0) (43) 

= J drdvr P(r){T,7r;t\0[ip(r + t)]\r,n;t) (44) 

= J drdit P(T){T,ir;t\0[-id n }\T,ir;t) (45) 

: (9{t)\d[<p]Mt)) , (46) 

6[p]=O[0\ , 0=-id^ , (47) 



where: 



in r, 7r-representation. In (|43|) - (j44|) the notation is symbolical, since the observable should be 
properly included in the functional integral defining the ket state, for example. 

Thus, a classical observable is represented by the corresponding function(al) of a suitably 
defined momentum operator. Furthermore, its expectation value at physical time t is repre- 
sented by the effective quantum mechanical expectation value of the corresponding operator 
with respect to the physical-time dependent state under consideration, which incorporates the 
weighted average over the proper times r, according to the distribution P. Not quite surpris- 
ingly, the evaluation of expectation values involves an integration over the whole r-parametrized 
"history" of the states. 

Furthermore, making use of the evolution operator U of Section 5, in order to refer ob- 
servables at different proper times t±, T2, . . . to a common reference point r, one can construct 
correlation functions of observables as well, similarly as in [S], for example. 

The most important observable for our present purposes is the classical Hamiltonian, H(tp), 
which enters the Hamiltonian constraint of a classical reparametrization invariant system. It 
is, by assumption, a constant of the classical motion. However, it is easy to see that also 
its quantum descendant, H((p) = H((p), is conserved, since it commutes with the effective 
Hamiltonian of ()37j) : 

— — d d 

[H,H] = H{-id^)n-u-—H{^)\ v= _ i9n -'K-uj-—H{^)\ v= _ idn H{-id^) 



—H( ( p)\ (p= _ ie ^u-—H( ( p)\ (p= ^ dn = , (48) 



due to the antisymmetric character of the symplectic matrix. Therefore, it suffices to implement 
the Hamiltonian constraint at an arbitrary time. 

Then, the constraint of the form Cu = H[(p] — e ~ may be incorporated into the definition 
of the states in ()24j) by including an extra factor 5(Ch) into the functional integral, and analo- 
gously for the adjoint states. Exponentiating the 5-function, we pull the exponential out of the 
functional integral, as before. Thus, we find the following operator representing the constraint: 

C = J d\ exp (zA(Sfa) - e)) = 5{C H ) , (49) 

which acts on states as a projector. Of course, a corresponding number of projectors should be 
included into the definition of the generating functional, see (I22J1 . for appropriate normalization 
of the states. 
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Supplementing (|42j ) - (j46|) by the insertion of the Hamiltonian constraint, the properly con- 
strained expection values of observables should be calculated according to: 

(0[<p];t) E = (9(t)\d[<p]d\9(ty , (50) 

which will deviate from the results of the previous definition. 

Finally, also the eigenvalue problem of stationary states, see (f4*Uj) - (|4*Tj) . should be studied 
in the projected subspace: 

HCW = EC\V) , (51) 
to which we shall return in the following examples. 



7 Examples of Emergent Quantum Systems 

7.1 Quantum Harmonic Oscillator from Classical One Beneath 

All integrable models can be presented as collections of independent harmonic oscillators. There- 
fore, we begin with the harmonic oscillator of unit mass and of frequency Q. The action is: 

^-/dt(^(^) 2 -^(g 2 -2e)) , (52) 

where A denotes the arbitrary lapse function, i.e. Lagrange multiplier for the Hamiltonian 
constraint, and e > is the parameter fixing the energy presented by this constraint. 

Introducing the proper time, r = / dt A, the Hamiltonian equations of motion and Hamil- 
tonian constraint for the oscillator are: 

d T q = p , d T p = -Vl 2 q , (53) 
i(p? + flV)-e = , (54) 

respectively. 

Comparing the general structure of the equations of motion Q15j) with the ones obtained 
here, we identify the effective Hamilton operator ([37)1 . while the constraint operator follows 
from flUIJ): 

H = -(vr 9 ^p - n 2 7r p (p q ) = -TT q (-id^ p ) + tt 2 7r p (-id nq ) , (55) 

C = 6(0 p 2 + fiV, 2 - 2e) = 6(d 2 p + tfd 2 q + 2e) , (56) 

respectively. Here we employ the convenient notation (p a = ((p q ;(p p ), and correspondingly 
7r a = (n^Tip), = (d nq ;d np ). Further simplifying this with the help of polar coordinates, 
7c q = —Qpcos<f) and tc p = psin0, we obtain: 

n = nl z = -iVLdj, , (57) 

C = 5(A 2 + 2e) = 6(d p 2 + p~ x d p + p~ 2 d 2 + 2e) , (58) 

where L z denotes the z-component of the usual angular momentum operator and A 2 the Lapla- 
cian in two dimensions. 
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We observe that the eigenfunctions of the eigenvalue problem posed here factorize into a 
radial and an angular part. The radial eigenf unction, a cylinder function, is important for 
the calculation of expectation values of certain operators and the overall normalization of the 
resulting wave functions. However, it does not influence the most interesting spectrum of the 
Hamiltonian. 

In the absence of the full angular momentum algebra, we discretize the angular derivative. 
Then, the energy eigenvalue problem consists in: 

H^n) = -i(nN/2rt)(ij(<j) n+1 ) - i){cj> n )) = £V(0n) , (59) 

with cf) n = 2nn/N , 1 < n < N, and the continuum limit will be considered momentarily. 
A complete orthonormal set of eigenfunctions and the eigenvalues are: 

<M0«) = N- 1 / 2 exp[t{m + 5)<f) n ] , l<m< A^ , (60) 
E m = i(nN/2n)(l -exp[2ni(m + 5)/N]) (61) 

Q(m + 5) , meN , (62) 

where 5 is an arbitrary real constant. 

Obviously, the freedom in choosing the constant 5, which arises from the regularization 
of the Hamilton operator, is very wellcome. Choosing 5 = —1/2, we arrive at the quantum 
harmonic oscillator, starting from the corresponding classical system. Thus, we recover in a 
straightforward way 't Hooft's result, derived from an equivalent cellular automaton j^J. See 
also [2] for the completion of a similar quantum model. In the following example we will 
encounter one more model of this kind and demonstrate its solution in detail. 

Here, and similarly in following examples, the eigenvalues are complex, with the real spec- 
trum only obtained in the continuum limit. This is due to the fact that we discretize first-order 
derivatives most simply, i.e. asymmetrically. It can be avoided easily by employing a symmetric 
discretization, if necessary. 

We find it interesting that our general Hamilton operator (|37jl does not allow for the direct 
addition of a constant energy term, while the regularization performed here does. 



7.2 Quantum System with Classical Relativistic Particle Beneath 

Introducing proper time as in Section 2, but leaving the extradimensions for now, the equations 
of motion and the Hamiltonian constraint of the reparametrization invariant kinematics of a 
classical relativistic particle of mass m are given by: 

d T q^ = m'V , d T p^ = , (63) 
p ■ p — m 2 = , (64) 

respectively. Here we have ip a = (q°, . . . ,q 3 ;p°, . . . ,p 3 ), a = 1, . . . , 8; four-vector products 
involve the Minkowski metric, g^ v = diag(l, —1,-1,-1). 

Proceeding as before, we identify the effective Hamilton operator: 

H = -m~ l Tt; q ■ (p p = -m~ l -K q ■ (~id np ) , (65) 
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corresponding to ()37|): the notation is as introduced after ()56)). however, involving four- vectors. 
Furthermore, the Hamiltonian constraint is represented by the operator: 

C = 5(0 2 - m 2 ) = 5{dl + m 2 ) , (66) 

following from (J49|) . 

After a Fourier transformation, which replaces the variable Tt q by a derivative (four-vector) 
+id x , and with 7r p = x, the Hamiltonian and constraint operators become: 

H = -m- x d x ■ d x , C = 5(d? + m 2 ) , (67) 

respectively. 

Similarly as in the harmonic oscillator case, the eigenvalue problem is properly defined by 
solved by discretizing the system on a hypercubic lattice of volume L 8 (lattice spacing I = L/N) 
with periodic boundary conditions, for example. Here we obtain the eigenfunctions: 

ipk x M( x n, x n ) = N' 1 exp[i(k x + 5 X ) ■ x n + i(k x + 5 X ) ■ x n ] , (68) 

with coordinates x% = In^ and momenta k% = 27r/c M /L, with 1 < n M , < N, and where 5% are 
arbitrary real constants, for all /i = 0, . . . , 3 (analogously x%, k£, 8%). 
The energy eigenvalues are: 

E kxM = -m-H- 2 ((expW(k x + 5 x ) }-l)(exp[il(k x + 5 x )°}-l) (69) 
- f>xp[z/(^ + 5 x y] - l)(exp[il(h + - 1)) 

= m -\k x + 8 x )-(k x + 5 x ) + 0(l) , (70) 

where is L is kept constant in the continuum limit, / — > 0. Furthermore, in this limit, one finds 
that the Hamiltonian constraint requires timelike "on-shell" vectors k x , obeying (k x +S x ) 2 = m 2 , 
while leaving k x unconstrained. 

Continuing, we perform also the infinite volume limit, L — > oo, which results in a continuous 
energy spectrum in (|70j) . We observe that no matter how we choose the constants 5 X ,5 X , the 
spectrum will not be positive definite. Thus, the emergent model is not acceptable, since it 
does not lead to a stable groundstate. 

However, let us proceed more carefully with the various limits involved and show that indeed 
a well-defined quantum model can be obtained. For simplicity, considering (1 + 1) -dimensional 
Minkowski space and anticipating the massless limit, we rewrite (fTOj) explicitly: 




where we suitably rescaled and renamed the constants and the momenta, which run in the range 
1 < k 1 , k 0,1 < N = 2s + 1. Furthermore, we incorporated the Hamiltonian (on-shell) constraint, 
such that only the positive root contributes: k° + 5° = \k l + 5 l \ + 0(m 2 ) = — (k 1 + 5 1 ) +0(m 2 ). 
This can be achieved by suitably choosing 5 0,1 . 

In fact, just as in the previous harmonic oscillator case, the choice of the constants is crucial 
in defining the quantum model. Here we set: 
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This results in the manifestly positive definite spectrum: 

E(s z , s /) = (^) 2 fe + s + I) + l) (>° + a + b + (si + s + ±) + l) + 0(m) , (73) 



mV ^ 2' 



with (half)integer quantum numbers s z , s^' 1 , all in the range — s < s z < s, replacing k l , k ' 1 . 

Recalling the algebra of the SU(2) generators, with S z \s z ) = s z \s z ) in particular, we are led 
to consider the generic operator: 

h = S z + s+ l - , (74) 

i.e., diagonal with respect to |s 2 )-states of the (half)integer representations determined by s. 
In terms of such operators, we obtain the regularized Hamiltonian corresponding to ()73|) : 

H = (-^-) 2 (l + h + h + In + h(h + hj) + 0(m) , (75) 



which will turn out to be equivalent to three harmonic oscillators, including a coupling term 
plus an additional contribution to the vacuum energy. 

A Hamiltonian of the type of h has been the starting point of 't Hooft's analysis which 
we adapt for our purposes in the following. 

Continuing with standard notation, we have S 2 = S 2 + S 2 + S 2 = s(s + 1), which suffices 
to obtain the following identity: 

h=-^—(S 2 + S 2 + - + h 2 ) . (76) 
2s + 1 v x v 4 ' v 1 

Furthermore, using S± = S x ± iS y , we define coordinate and conjugate momentum operators: 

q=±(aS- + a*S+) , p + , (77) 

where a and b are complex coefficients. Calculating the basic commutator with the help of 
[<S+, S-] = 2S Z and using (|7%J). we obtain: 

provided we set 5s (a* b) = —2/ (2s + 1). Incorporating this, we calculate: 

S 2 + S y 2 = (2S ± 1)2 (\a\ 2 f + \b\ 2 q 2 - (5sa -Sb + fta- W>){q,p}) . (79) 
In order to obtain a reasonable Hamiltonian in the continuum limit, we set: 

a = i , , b= , fi = (^-) 2 . (80) 

1 8 + 1/2 Js + 1/2 V™^ 



Then, the previous (fTE|) becomes: 



™ = ^ + ^¥ + ^ Tp (itf + (^) , (81) 
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reveiling a nonlinearly modified harmonic oscillator Hamiltonian, similarly as in [21 Ej. 

Now it is safe to consider the continuum limit, 2s + 1 = iV — > oo, keeping y/mL and Q 
finite. This produces the usual g, j3-commutator in (|78jl for states with limited energy and the 
standard harmonic oscillator Hamiltonian in (|81j) . 

Using these results in (|75jl . the Hamilton operator of the emergent quantum model is ob- 
tained: 

n = n + 1 E fc 2 + ^/) + ^#i 2 + ^i 2 ) E fe 2 + ^/) , (82) 

i=i,o,i j=o,i 

where the massless limit together with the infinite volume limit is carried out, m — > 0, L — > oo, 
in such a way that Q remains finite. 

The resulting Hamiltonian here is well defined in terms of continuous operators q and p, 
as usual, and has a positive definite spectrum. The coupling term might appear slightly less 
unfamiliar, if the oscillator algebra is realized in terms of bosonic creation and annihilation 
operators. 

We previously calculated the matrix elements of operators q,p with respect to the SU(2) 
basis of primordial states in a similar case, showing that localization of the quantum oscillator 
has little to do with localization in the classical model beneath |IJ |2] . 

Finally, we remark that had we chosen 5 0,1 = 5 0,1 = 1/2, instead of (|72jl. then a relative 
sign between terms would remain, originating from the Minkowski metric, and this would yield 
the Hamiltonian 7i oc (1 + h)(h — hi), which is not positive definite. Similarly, any symmetric 
choice, <f 0,1 = 5 0,1 = 5 would suffer from this problem. 

This raises the important issue of the role of canonical transformations, and of symmetries 
in particular. It is conceivable that symmetries will play a role in restricting the present 
arbitrariness of the regularization defining a quantum model. We will address further aspects 
of this in the following section. 



8 Remarks on (Non)Integrable Interactions 

We resume our discussion of general features of the emergent quantum mechanics. Specifically, 
let us consider a classical system with n degrees of freedom, for example, a chain of particles 
with harmonic coupling and anharmonic potentials. Denoting the phase space variables by 
if a = (Q,P), where Q,P are n-component vectors, we assume for definiteness a Hamiltonian 
of the form, H(ip) = (1/2)P 2 + V(Q), i.e. with a kinetic term which is simply quadratic in the 
momenta. 

In this case, following (pfTJl and (J49|) . and with: 

Q = -i&KQ , P = -id^p , (83) 
the Hamiltonian and constraint operators, respectively, are given by: 

H = -n Q -P + n P -V'(Q) , (84) 
C = 6(±P 2 + V(Q)-e) , (85) 

where, of course, V'(Q) = VqV(Q), and the wave function is considered as a function ip(irp, tcq) 
of the indicated vectors. 
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The previous oscillator example suggests to perform a Fourier transformation to variables 
x, y, such that the eigenvalue problem becomes: 

HiJ>(x,y) = (x ■ (-id y ) - V'(y) ■ (-id x ))i/;(x,y) = E^{x,y) , (86) 

while the constraint operator equation turns into an algebraic constraint: 

Cif>(x, y) = 5(^x 2 + V(y) - e)i/>{x, y) = . (87) 

In agreement with the general result ()48|) . we easily confirm here that TiCip = CTiip. The 
constraint equation then simply states that the phase space variables are constrained to a 
constant energy surface of the underlying classical system. 

The first order quasi-linear partial differential equation (J86)) can be studied by the method 
of characteristics [201 ■ Thus, one finds one equation taking care of the inhomogeneity (right- 
hand side), which can be trivially integrated. Furthermore, the remaining 2n equations for the 
characteristics present nothing but the classical Hamiltonian equations of motion. 

It follows that integrable classical models can (in principle) be decoupled in this context of 
the characteristic equations by canonical transformations. This assumes that we can apply them 
freely at the pre-quantum level, which might not be the case. It would lead us essentially to the 
collection of harmonic oscillators mentioned at the beginning of Section 7.1, and corresponding 
quantum harmonic oscillators as studied there. 

Classical crystal-like models with only harmonic forces, or free field theories, respectively, 
will thus give rise to corresponding free quantum mechanical systems here. These are con- 
structed in a different way in Presumably, the (fixing of a large class of) gauge transforma- 
tions invoked there can be related to the existence of integrals of motion implied by integrability 
here. In any case, we conclude that in the present framework truly interacting quantum (field) 
theories might be connected with nonintegrable deterministic systems beneath. 

Furthermore, we emphasize that the Hamiltonian equations of motion preclude motion 
into classically forbidden regions of the underlying system. Nevertheless, quantum mechanical 
tunneling is an intrinsic property of the quantum oscillator models that we obtained, as well as 
of the anharmonic oscillator example considered in the following. Similarly, spreading of wave 
packets is to be expected in the latter case. 

In order to demonstrate additional features of the eigenvalue problem of (|86|). we concentrate 
on one degree of freedom with phase space coordinates p, q and with a generic anharmonic 
potential. 

Since the potential depends only on q, by locally stretching or squeezing the coordinate, 
i.e. by an "oscillator transformation" q = f(q), we can bring it into oscillator form, such that 
V(f(q)) = (l/2)q 2 . Implementing this type of transformation, the equations for one degree of 
freedom are: 

(pdg - qd p )^j(p, q) = Eip(p, q) , (88) 

S(l(p 2 + q 2 )-e)^(p,q) = , (89) 

with /' denoting the derivative of /. This is very much oscillator-like indeed and, once more 
employing polar coordinates, we obtain: 



(90) 



20 



5(p 2 -e)^(p,0) = O . (91) 

The eigenvalue problem seems underdetermined. As it stands, it would give rise to an unbound 
continuous spectrum, with no groundstate in particular. 

This apparent defect persists for any number of degrees of freedom. However, as we have 
seen already, sense can be made of the Hamilton operator by a suitable regularization, especially 
by discretizing the phase space coordinates. The principles of such regularization we still do 
not know, other than either preserving or intentionally breaking symmetries. 

8.1 An Anharmonic Oscillator 

It is worth while to consider one more example, a onedimensional system with Hamiltonian: 

H = ^p 2 + V \q\ , (92) 

in order to demonstrate the subtleties associated with regularization. For the linear poten- 
tial, the coordinate dependence of the operator on the left-hand side of (J9T?|) is mild, since 
/'(psin0) = V^~ 1 psin0, and the eigenvalues of its discretized counterpart can be found as 
follows. 

Conveniently discretizing the angular variable as <p n = (27rn/N) + (37r/2), 1 < n < N, the 
eigenvalue equation becomes: 

N 

H (l - Acos(27m/iV)) = 1 , A = 2niy/eE/NV . (93) 

n=l 

Setting A = 2z/(l + z 2 ), and employing a known identity for the finite product arising here, 
the eigenvalue equation can be transformed into: (1 + z N ) 2 = (1 + z 2 ) N . With hindsight, we 
choose iV = 2(4N' + 1) and set z = +\/u — 1, to obtain: 

u N/2 + {1 _ u) N/2 = 1 (94) 

This equation has the nice property that, if u is a solution, then so is 1/u. The location of the 
solutions u — 0, l,oo suggests to look for further solutions in the form of u = exp2ia. Thus, 
combining the equations for u and 1/u, we arrive at the transcendental equation: 

2sin(aiV/2) = (2sina) 7V/2 . (95) 

From the multitude of its solutions, due to periodicity, we need to find iV solutions of (|9*3^) . 

Closer inspection shows that, in the limit of large N, solutions of (J95j) consist essentially 
of those zeros of sin(aiV/2) which lie inside the intervals [mr — 7r/6,mT + vr/6], with integer n. 
Thus, positive energy solutions will be obtained momentarily from: 

E± = NV exp(i^ ± ~ia) ((2/e) sin(±«)) 1/2 , (96) 

where either "+" or "— " has to be chosen consistently, corresponding to the solutions coming 
in pairs exp ±2ia. 

A remark is in order here. Considering only the positive root above, z = +a/m — 1, we 
avoided negative energy solutions. However, there is a price to pay: careful counting reveils 
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that the positive energy spectrum is doubly degenerate. The finite positive part is obtained 
from (jUnj), incorporating a = 2iim/N: 

1 /2 

E = NV exp(3i7rm/N)((2/e) sin(27rm/iV)) , < m < m < N/12 , (97) 

y (m + m -l) 1/2 , m G N , (98) 

where m is an arbitrary constant, within the allowed range, which defines the zeropoint en- 
ergy of the emergent quantum model. The continuum limit is to be taken such that Vo = 
Vo(47rA r /e) 1 / 2 stays finite. The additional solutions can be chosen in a way that their real parts 
move to +oo, as N — > oo. 4 

We remark that the spectrum of (}9~H|) differs from the one obtained for the same potential 
in standard quantum mechanics, where WKB yields: E oc (m — l/4) 2//3 . 

Summarizing, the various illustrated features promise to make genuinely interacting models 
quite difficult to analyze. We hope that more interesting results will be obtained with the help 
of spectrum generating algebras or some to-be-developed perturbative methods. 



9 Conclusions 

We pursue the view that quantum mechanics is an emergent description of nature, which 
possibly can be based on classical, pre-quantum concepts. 

Our approach is motivated by a construction of a reparametrization-invariant time. In turn, 
this is based on the observation that "time passes" when there is an observable change, which 
is localized with the observer. More precisely, necessary are incidents, i.e. observable unit 
changes, which are recorded, and from which invariant quantities characterizing the change of 
the evolving system can be derived. 

We recall the model of [2], invoking compactified extradimensions in which a particle moves 
in addition to its relativistic motion in Minkowski space. We employ a window to these extradi- 
mensions, i.e., we consider a quasi-local detector which registers the particle trajectory passing 
by. Counting such incidents, we construct an invariant measure of time. 

A basic ingredient is the assumption of ergodicity, such that the system explores dynamically 
the whole allowed energy surface in phase space. This assures that there are sufficiently frequent 
observable incidents. They reflect properties of the dynamics with respect to (subsets of) 
Poincare sections. Roughly, the passing time corresponds to the observable change there. 
Then, the particle's proper time is linearly related to the physical time, however, subject to 
stochastic fluctuations. 

Thus, the reparametrization-invariant time based on quasi-local observables naturally in- 
duces stochastic features in the behavior of the external relativistic particle motion. Due to 
quasi-periodicity (or, generally, more strongly irregular features) of the emerging discrete time, 
the remaining predictable aspects appear as in unitary quantum mechanical evolution. 

In reparametrization-invariant, "timeless" single-particle systems, this idea has been realized 
in various forms PUI2]- Presently, this has led us to assume the relation between the constructed 
physical time t and standard proper time r of the evolving system in the form of a statistical 
distribution, P(r; t) = P(r — t), cf. (J20J). Here we assume that the distribution is not explicitly 

4 Corresponding eigenfunctions, i.e. iV-component discrete eigenvectors, are obtained by evaluating products 
of the kind appearing in (|93[) . with k — 1 < N — 1 factors for the kth component and a constant for k = 1. 
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time-dependent, which means, the physical clock is decoupled from the system under study. 
We explore the consequences of this situation for the description of the system. 

We have shown how to introduce "states", eventually building up a Hilbert space, in terms 
of certain functional integrals, (|24 p -([25| h which arise from the study of a suitable classical 
generating functional. The latter was introduced earlier in a different context, studying classical 
mechanics in functional form We employ this as a convenient tool, and modify it, in 

order to describe the observables of reparametrization-invariant systems with discrete time 
(Section 6). Studying the evolution of the states in general (Section 5), we are led to the 
Schrodinger equation, (|39j). However, the Hamilton operator (|37|) has a non-standard first- 
order form with respect to phase space coordinates. 

The choice of boundary conditions of the classical paths contributing to the generating 
functional plays a crucial role and deserves better understanding. 

Furthermore, illustrating the emergent quantum models in various examples, we demon- 
strate that proper regularization of the continuum Hamilton operator is indispensable, in order 
that well-defined quantum mechanical systems emerge, with bounded spectra and a stable 
groundstate, in particular. Most desirable is a deeper undertanding of this mapping between 
the continuum Hamilton operator, which is straightforward to write down, given a classical 
pre-quantum system, and the effective quantum mechanics, which emerges after proper regu- 
larization only. Especially, limitations imposed by symmetries and consistency of the procedure 
need further study. 

It is a common experience that the preservation of continuum symmetries through dis- 
cretization is difficult, for example, see [T^l [T3J [THl H7| . We wonder, whether other regularization 
schemes are conceivable. The possiblity, mentioned after (|39|). that the need for regularization 
is an artefact of decoupled clock degrees of freedom deserves further study. 

We find that truly interacting quantum (field) theories might be connnected to nonintegrable 
classical models beneath, since otherwise the degrees of freedom represented in the stationary 
Schrodinger equation here, can principally be decoupled by employing classical canonical trans- 
formations. 

Finally, we come back to the probabilistic relation between physical time and the evolution 
parameter figuring in the parameterized classical equations of motion, which is the underlying 
raison d'etre of the presented stroboscopic quantization. One would like to include the clock 
degrees of freedom consistently into the dynamics, in order to address the closed Universe. 
This can be achieved by introducing suitable projectors into the generating functional. Their 
task is to replace a simple quasi-local detector which responds to a particle trajectory passing 
through in Yes/No fashion; by counting such incidents, an invariant measure of time has been 
obtained before [HI2|- I n a more general setting, this detector /projector has to be defined in 
terms of observables of the closed system. In this way, typical conditional probabilities can be 
handled, such as describing 'What is the probability of observable X having a value in a range 
x to x + 5x, when observable Y has value yT' . Criteria for selecting the to-be-clock degrees 
of freedom are still unknown, other than simplicity. Most likely the resulting description of 
evolution and implicit notion of physical time will correspond to our distribution P(r; t) of 
(|2U|) . however, now evolving explicitly with the system. We leave this for future study. 

The stroboscopic quantization emerging from underlying classical dynamics may be ques- 
tioned in many respects. It might violate one or the other assumption of existing no-go theorems 
relating to hidden variables theories. However, we believe it is interesting to learn more about 
working examples, before discussing this. Unitary evolution, tunneling effects, and spreading 
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of wave packets are recovered in this framework. Interacting theories remain to be explored. 
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